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TECHNICAL MEMORANDUM NO. 1100 

INFINITESIMAL CONICAL SUPERSONIC FLOW* 

By Adolf Busem&hn 


Conical flow fields .- Real flows always occur in 
three-dimensional space. In calaulatiiig a flow, however, 
one will greatly appreciate It If there are only two 
essential coordinates to deal with. Flows of this .kind', 
limited to two coordinates, form the plane flow and the 
flow of axial symmetry. The space which Is filled out 
by the streamlines Is represented in planes parallel to 
these lines; they contain certain streamlines to their 
whole extent. In conical flow fields, however, the 
streamlines are cut through slantingly so that each 
streamline Is contained in the plane but appears there 
as a point only. These relations are made clear In fig- 
ure 1. If the friction Is neglected the shape of the 
body leads one to expect a pattern that can be increased 
or decreased geometrically. The fixed point P and the 
direction of the tliree spatial axes, X, y, and & remain 
the same. All essential characteristics of the flow and 
the shape of the body can be inferred from the plane 
z = 1 . a pi. mo z = 2 \/ould, if lictuucuS were doubled,, 
show identical values for gas conditions and velocities. 
The Isobar planes In the space x, y, 2 are of conical 
shape and have the cone vertex P; therefore these flows 
shall be called abbrev^atedly conical flow fields. 

Infinitesimal differences in pressure .- In figure 1 
there shall be one more limitation r or the general conioal 
flow field, namely that the body disturbs the parallel 
flow only to a slight degree. So the conical isobars 
reflect over- and under-pressures differing Infinitesi- 
mally from the pressure of the parallel flow. This two- 
fold limitation, to conical and infinitesimal, is not 
actually very stringent Insofar as In the class of poten- 
tial flows there are present only the conical fields of 
axial symmetry and the Infinitesimal conical fields. All 
other conical flows are affected by rotation. The Infini- 
tesimal supersonic flows, however, also excel In another 
way: the superposition of fields with different fixed 

? oints P is permitted in spite of the fact that the inf- 
erential equations ordinarily are not linear; thus the 
applicability Is broadened most gratifylngly. 

*" Infinitesimal© kogeTIge tibersohallstrSmung." 
Deutschen Akademle der Luftf ahrtforshung, I 9 I 42 -I 4 . 5 , p. Ij. 55* 
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Differential equation.- There are two ways to limit 
the dirrerential equation for the potential ir^ conical 
fields to small additional velocities u, v, w, and to 
limit the differential equation for nearly parallel 
spatial flows to conical fields; the first is the his- 
torical one, the second, however, the simpler one. 
Therefore here the second one is chosen. As is well 
known, the linearized differential equation in the 
space x, y, z, for the additional potential cp over a 
basic velocity W in the direction of" the axis z reads 
if the gas has the sonic velocity a: 



The coordinates £ and tj of the conical current cor- 
respond to the spatial coordinates' x and y in the 
plane 2=1; 

£ = £ and n = \ (2) 

Z 3 

The additional potential op increases on each ray through 
the fixed point P in proportion to the distanoe. There- 
fore, the potential divided by z is invariant on a 
single ray, and provides the potential of the conical 
flow; 


v = — <P (x,y,z )_^ (?) 

The additional velocities u, v, w arc the derivatives 
of the formor potential. 


u - <?x = 

\7 = <p z = % - l 


v = "V = \ 

* Pg “ T] * Xtj 






J 


(b) 


The differential equation for the new potential y is 
determined from the old differential equation, and one 
obtains : . 
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with 





( 5 ) 


It certainly Is gratifying to recognize In the type of 
this differential equation an old acquaintance from the 
plane gas flow: for thostroam function of the plane flow 
transformed according to Legendre and superimposed over 
the components <9f the current density produces exactly 
the same differential equation. In ordinary gases, how- 
ever, the denominator A 3 is a local function; but there 
Is a special gas with rectilinear adiabatics In the 
pre sour e -volume -diagram In which, as required, this 
denominator also remains constant. This gas is a special 
favorite where it is a mere question of numerical 
calculations . 


Regions of influence .- The spatial differential 
equation of the gas ^low at supersonic velocity is of 
hyperbolic character, as shown in equation (1). That 
means: each point of the flow dominates a conical range 

opening downstream; each locus, on the other hand, is 
dominated solely by those points which are situated in 
the cone prolonged backward and opening upstream. Here- 
with the relationships are divided definitely among the 
three possibilities • superior, subordinate, and independent. 
Macli*s cones in the supersonic flow considered as regions 
of disturbance of a small trial body make this fully com- 
prehensible in a physical sense. It must seem odd at 
first that the dependencies of the general spatial flow 
are widened as soon as one proceeds to a more limited 
spatial flow. But the above-mentioned differential equa- 
tion shows thDt Inside of the circle with the radius A 
there prevails the elliptic character. 


This behavior is easily explained by the fact that 
all points of a ray starting from P are comprehended 
as a whole. The relation of dependencies of two rays 
results from the dependencies of the single points. Only 
the char actsrls tic "independent" appears uniformly In 
certain cases for all pairs of points (P Itself is 
excluded). The combination superior and independent 
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becomes superior; subordinate, and independent become 
subordinate. But if there are pairs of points of all 
kinds on the rays, then the rays are subject to the new 
characteristic "reciprocally dependent." Rays of this 
kind fill out the interior of Mach's cone starting from 
point P. 

Characteristics.- Mach's cone starting from P 
intersects the plane z = 1 on the circle having the 
radius A. In the field outside of this' cone, i. e., 
outside of the circle In the Intersecting plane, one 
gets rectilinear characteristics of the differential 
equation (5) which are tangents of the circle. In fig- 
ure 2 this Is demonstrated by two wires, a and b. 

The wire b Is bent slightly upstream in order not to 
exclude cases of this kind. The range of disturbance 
results from the sum of all of Mach's cones starting from 
all points of the wire. It is immed? ately obvious that 
only the circle with the radius A and its tangents can 
form the boundaries of the area of disturbance. Outside 
of Llach's cone starting from point P these character- 
istics settle all questions; they can be traced back to 
the plane case with a transverse component of the velocity. 
The essential and different part of the conical fields, 
therefore. Is concerned with the convex surface of Mach's 
cone starting from point P, and with Its Interior. 

Tschapligin* s illustration.- In the plane of Inter- 
section z =~1 we find iris l ae of the circle with the 
radius A the elliptic charactor of the differential 
equation (5). Near the center the differential equation 
of the- potential theory is valid; In plane cases, this 
equation can be satisfied by analytic functions of the 
complex variable.- In this circle, therefore, there only 
exists a mutual dependence but not yet a full equivalence 
of all loci. This is not surprising, because the analytic 
continuation of the plane reaches to the outer range of 
the circle. Tschapligin, however, lia3 devised a geomet- 
rical construction which so distorts the field inside the 
circle that equivalence regarding the differential equa- 
tion will result. As figure 3 shov/a this distortion is 
attained by transferring the plane z =1, with the 
complex variable £ = £ + Ip,, through parallel projec- 
tion to a sphere with the radius A, and by then pro- 
jecting it from a pole of th9 sphere on to a plane In 
the distance. One will easily recognize that only the 
interior of the circle with tho radius A will be . 
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depicted; first it will be delineated from the lower 
half-sphere on the. interior of the unit circle of the 
new plane with the new complex variable - e; a second 
time it will go from the upper half-sphere on to the 
outer field of the unit circle. In these coordinates 
one can use analytic functions for the solutions. 


SOLUTION OP THE DIFFERENTIAL EQUATION 


For each of the velocity components u, v, and w 
one can equate the real part of an analytic function f(e). 
It will serve the purpose best to set up the equation 
for the component w, because then the more closely 
related components u and v can be calculated Jointly* 

w = A • Re (f ( e ) ) or w + Is = A • f (e) (6) 


The completion represented here by s is for the time 
being completely meaningless. According to Tschapligin 
there then results the complex velocity: 


co = u + iv 



( 7 )- 


The pressure in the current, with the aid of the den- 
sity p, results from the velocity components as follows 



(8) 


The right function f(e) Is to be selected with the aid 
of the boundary conditions. 

Boundary conditions .- The outside of Mach's cone Is 
superior to Mach's cone Itself. Therefore, first, those 
velocities u, v, and w on the circle of the £ -plane 
(and therefore on the uniform circle of the e-plane) 
thab result from the outer field must be ascertained. 



6 


NAG A TM No. 1100 


If the body does not protrude anywhere out of Mach’s 
cone, the values u = v = w= 0 on the uniform circle 
are given. If on the contrary no part of the- body Is 
Inside of Mach’s cone, the values of w are to be 
represented by an analytic function f ( e ) free of 
singularities with the given boundary values on the 
circle. If f ( c ) does not produce a stationary value 
df = 0, then u and v according to equation (7) will 
have a logarithmic singularity at zero. The many-leaved 
function can be selected in a ■unique way by using radial 
intersections with the boundary values of u and v on 
the uniform circle. The radial intersections produce 
rotational layers, as is physically to be expected from 
a lifting surface. 

Impermeable boundaries of the body can be trans- 
ferred into the £ -plane at the same time. They mu3t be 
streamlines in the field of the relative velocity: 


fi rol “ “ + fi2a) - 2AiV£ (9) 

This condition is not always easy to comply with. How- 
ever, if the body possesses rectilinear surface elements 
passing near zero, the otherwise meaningless imaginary 
part s of the function f(e) will remain constant on 
these elements. If bhe straight part goes ovor zero, a 
stationary value for f is to be stipulated at zero. 
Conditions of this kind arc especially agreeable. Prom 
the pressure equation (8) conditions applicable to cases 
of given pressures or of given lifts are to be understood. 

The disappearance of the real or of the imaginary 
part of f on certain lines because of symmetry can be 
attained in the woll known way by reflexion, as the 
examples will show. 


Examples 

1, The circular cone In the straight flow 

Pop the only axial- symmetrical case, 1. e., the 
circular cone with art infinitesimal apex angle, the 
right solution Is, of course, given by the statement 

* ine 


w + 13 = C 
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The pressure on the convex surface of the cone results 
In the known way (fig* I4.) and .conforms, with v, Karman*s 
values and mine. 

2. The circular cone in oblique flow 

One succeeds, with the aid of the relative velocity 
according to equation (9)* in solving the circular cone 
in oblique flow. Herein the apex angle and the angle of 
attack may, though infinitesimal, yet bear a relation to 
each other. The solution la shown in figure 5« If °ne 
makes the angle of obliquity y zero, one gets again the 
oircular cone in straight flow. If one makes the apex 
angle 2(3 disappear, one gets the pressure distribution, 
of a circular cone in an incompressible current. The 
comparison with Ferrari IS rendered somewhat difficult 
by the fact that Ferrari measures the velocity field per- 
pendicular to the cone axis while it is here perpendicular 
to the wind direction. If the system of coordinates Is 
rotated adequately, the conformity is complete. 

3. Tip of a rectangular plato 

If a plane rectangular plate of infinitesimal thick- 
ness is placed in a flow perpendicular to the front edge 
with an inf initeoimal angle of attack, and if the velocity 
field is needed only up to the roar edge of the plate, 
one can place the fixed point P at the right corner 
point of the front edge. On the supposition of an 
Infinitesimal angle of attack y (with the x-axis 
forming the axis of rotation) the pressure distribution 
will be represented on the quarter plane between the 
positive axis (z) and the' negative axis (x). For the N 
plane z = 1 tho section of the body, except for infini- 
tesimal distances, is then rendered by the negative real 
axis. Let the reduced pressure above the plate and the 
increased pressure below the plate be adjusted to a unit 
value outside of Mach»s cone. These values hold on the 
boundary circle. On the left half of the unit circle of 
the £ -plane corresponding values for w are then to be 
assigned. On the right semicircle the outer field is 
undisturbed; h3re w = 0. For reasons of symmetry the 
value w = 0 must also result on tho positive, real axis. 
Along the negative ro$l axis, on the contrary, 
s = lm(f(e)) must bo fixed because of the fixed radial 
boundary. Since 3 Is given only up to one constant, 
one can demand her9 s = 0. All conditions can be 
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attained by reflexion If one undertakes a preliminary 
conformal mapping on the plane U = -\/ 7 . 

The solution is represented in figure 7* Figure 8 
shows the pressure distribution on both edges of a rec- 
tangular plate. 


4. Supporting triangle 

Every tv/o radii starting from P form a triangular 
plane as far as the plane z = 1, when all points are 
connected. Because of the required inf initesimal dis- 
turbance of the parallel current, hov/ever, the angle of 
attack must be infinitesimal, so that the plane of the 
two rays will nearly pass through the z-axis. Such 
triangles are possible completely inside of Mach's cone, 
completely outside, and mil- and bi-laterally protruding. 
Here we shall only consider the simplest case of the 
supporting triangle outside of loach's cone, although all 
other cases can be easily integrated. 

Figure 9 shows this supporting triangle. The 
velocity component w which predominantly influences 
the pressure is different from zero only on the short 
arcs between £ x and £ 3 as also £ 4 and £ a . The 
value zero results from the undisturbed state on the 
right, and also on the left because of the pressure 
adjustment behind the triangle, when consideration is 
given to the symmetry with a positive and wath a nega- 
tive angle of attack. Figure 10 shows the relations in 
the J-plane. If one Intends to let the rear edge of 
the triangle travel while the front edge lies fixed, one 
will at first transfer only the points £ x and £ a into 
the e-plane. With suitable regulation there must result 
an increase of w from 0 to +ir at £3., and from 

-tt to 0 at 5 3 (If one moves on the circle In the 

direction of increasing angles). One can treat this 
part of the solution independently if one assumes a 
further singularity at zero. Physically speaking, one 
then has a uniformly loaded triangle between the front 
edge at and the z-axis. Inside of llach's cone, 

however > t is triangle is not flat, but is twisted to 
uniform or load,.- As soon as on<* superimposes at 
the rear edge a negatively loaded triangle and its Influ- 
ence between £ KS and the z-axis, the part behind the 

rear edge will no longer be supporting, and the singu- 

larity in w at the poin 1 ; zero of the e-plane disappears. 
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However, a vortex layer In the field u, v Is left. 
The partial solution in the e-plane is represented In 
figure 'll; ■ - 


5. Superposition of two conical flows 

The infinitesimal conical flows can bo superposod 
without having the fixed point in common as in figure i|. 
Therefore the relations In the plane plate can also be 
represented when the plate has more depth. Figure 12 
shows the isobars of the edge of the plate and also their 
superposition after Mach«s cones have overlapped. Dif- 
ferent boundary conditions for the partial solutions 
need be considered only when the cones roach the other 
edge of the plate. The disappearance of pressure along 
a straight line in just the distance at which the cones 
arrive at the other edge of the plate is remarkable. 

The positively loaded part of the plate ends here. Fig- 
ure 15 3hows the lift distribution of the positively 
loaded part In perspective representation. 

To find the volocity field behind a rectangular 
plate of finite depth oac can annul the supporting pres- 
sure differences of a plate of Infinite depth by conical 
fields having the apices on the rear edge of the plate. 

If one superposes a negatively loaded plate shifted 
inf Initosimally in the direction of the z-axi3 one 
obtains a supporting line as a limiting case of the 
supporting strip. The cases calculated by Schlichting 
according to FrandtHs method nre obtained in this way. 
Hei’e, too, the conformity is perfect, except for an 
error of sign in the calculation of the Integral equation. 


SUMMARY 


The calculation of Infinitesimal conical supersonic 
flows has been applied first to the simplest examples 
that havo also boon calculated In another way. Except 
for the discovery of a miscalculation In an older report 
there was found the expected conformity. The new method 
of calculation Is limited more definitely to tho 'conical 
case* but, as a compensation. It is much more convenient 
because the solution is obtained by analytic functions. 
The fundamental recognition that there the hyperbolic 
character is roplaced by the elliptic ono will lead to 
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more thorough investigation of conical fields as special 
cases in supersonic flows. Of course, one will be 
tempted to call the elliptic character seemingly elliptic 
only | for if one notches an indentation into a cone, the 
real hyperbolic character of the field of the flow down- 
stream from this indentation immediately becomes obvious 
again. However, if one traces the flow to a point very 
far behind the indentation there will appear inside of 


Mach * s cono a reciprocal relation between every two rays 
which will gradually restore the conical course of the 


flow. Instead of trying to produce the final flow by an 


infinite succession of hyperbolic dependencies it will 
be more expedient to consider special ' elliptic singu- 
larities at the ocints of disturbance. In those rela- 


tions I see the significance of 


conical field; the 


infinitesimal case represents only a first approximation 


to it . 


Translation by Mary L. Mahler 
and Robert T. Jones, 

National Advisory Committee 
for Aeronautics . 




Figs. 3,4 
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/ = Cine 

p = z * 2 

A = tgCc 



Figure 4. Circular cone in axial flow. 









